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Based on lectures given at the 2008 Theoretical Advanced Study Institute (TASI), 
I review here some aspects of the phenomenology of particle dark matter, including 
the process of thermal freeze-out in the early universe, and the direct and indirect 
detection of WIMPs. I also describe some of the most popular particle candidates 
for dark matter and summarize the current status of the quest to discover dark 
matter's particle identity. 



1. Evidence For Dark Matter 

A wide variety of evidence has accumulated in support of dark matter's 
existence. At galactic and sub-galactic scales, this evidence includes galac- 
tic rotation curves [1], the weak gravitational lensing of distant galaxies 
by foreground structure [2], and the weak modulation of strong lensing 
around individual massive elliptical galaxies [3] . Furthermore, velocity dis- 
persions of stars in some dwarf galaxies imply that they contain as much 
as ~ 10^ times more mass than can be attributed to their luminosity. On 
the scale of galaxy clusters, observations (of radial velocities, weak lens- 
ing, and X-ray emission) indicate a total cosmological matter density of 
Qm ~ 0.2 — 0.3 [4], which is much larger than the corresponding density in 
baryons. In fact, it was measurements of velocity dispersions in the Coma 
cluster which led Fritz Zwicky to claim for the first time in 1933 that large 
quantities of non-luminous matter are required to be present [5]. On cos- 
mological scales, observations of the anisotropics in the cosmic microwave 
background have lead to a determination of the total matter density of 
Om^^ = 0.1326 ± 0.0063, where h is the Hubble parameter in units of 
100 km/sec per Mpc (this improves to f^M^^ = 0.13581q;qo36 distance 
measurements from baryon acoustic oscillations and type la supernovae are 
included) [6]. In contrast, this information combined with measurements of 
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the light chemical element abundances leads to an estimate of the baryonic 
density given by f^s/i^ = 0.02273±0.00062 (Jls^^ = m267t°amall ifBAO 
and SN are included) [7, 6]. Taken together, these observations strongly 
lead us to the conclusion that 80-85% of the matter in the universe (by 
mass) consists of non-luminous and non-baryonic material. 

The process of the formation of large scale structure through the gravita- 
tional clustering of coUisionless dark matter particles can be studied using 
N-body simulations. When the observed structure in our universe [8] is 
compared to the resiilts of cold (non-relativistic at time; of structure for- 
mation) dark matter simulations, good agreement has been found. The 
large scale structure predicted for hot dark matter, in contrast, is in strong 
disagreement with observations. 

Although there are many pieces of evidence in favor of dark matter, it is 
worth noting that they each infer dark matter's presence uniquely through 
its gravitational influence. In other words, we currently have no conclusive 
evidence for dark matter's electroweak or other non-gravitational interac- 
tions. Given this, it is natural to contemplate whether, rather than being 
indications of dark matter's existence, these observations might instead be 
revealing departures from the laws of gravity as described by general rela- 
tivity. 

Since first proposed by Milgrom in 1983 [9], efforts have been made to 
explain the observed galactic rotation curves without dark matter within 
the context of a phenomenological model known as modified Newtonian 
dynamics, or MOND. The basic idea of MOND is that Newton's second 
law, F = ma, is modified to F = ma x /L((a), where is very closely 
approximated by unity except in the case of very small accelerations, for 
which /i behaves as fi = a/aQ. Applying the modified form of Newton's 
second law to the gravitational force acting on a star outside of a galaxy of 
mass M leads us to 

r = J — = ma/j., (1) 



which in the low acceleration limit (large r, a<^ ao) yields 

y/GMao 



(2) 



r 

Equating this with the centrifugal acceleration associated with a circular 
orbit, we arrive at 

^^ = ^ ^ v = {GMao)'/'. (3) 
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In other words, MOND yields the prediction that galactic rotation 

curves should bcconic flat (independent of r) for sufficiently large orbits. 
This result is in good agreement with galaxy-scale observations for a value 
of oo ~ 1.2 X 10~^° m/s^, even without the introduction of dark matter. 
For this value of ao, the effects of MOND are imperceptible in laboratory 
or Solar System scale experiments. 

MOND is not as successful in explaining the other evidence for dark 
matter, however. In particular, MOND fails to successfully describe the 
observed features of galaxy clusters. Other evidence, such as the cosmic mi- 
crowave background anisotropics and large scale structure, are not generally 
able to be addressed by MOND, as MOND represents a phenomenological 
modification of Newtonian dynamics and thus is not applicable to ques- 
tions addressed by general relativity, such as the expansion history of the 
universe. Efforts to develop a viable, relativistically covariant theory which 
yields the behavior of MOND in the non-relativistic, weak- field limit have 
mostly been unsuccessful. A notable exception to this is Tensor- Vector- 
Scalar gravity, or TeVeS [10]. TeVeS, however, fails to explain cluster-scale 
observations without the introduction of dark matter [11]. This problem 
has been further exacerbated by recent observations of two merging clusters, 
known collectively as the bullet cluster. In the bullet cluster, the locations 
of the baryonic material and gravitational potential (as determined using 
X-ray observations and weak lensing, respectively) are clearly spatially sep- 
arated, strongly favoring the dark matter hypothesis over modifications of 
general relativity [12]. 

2. The Production of Deirk Matter in the Early Universe 

The nucleons, electrons and neutrinos that inhabit our universe can each 
trace their origin back to the first fraction of a second following the Big 
Bang. Although we do not know for certain how the dark matter came 
to be formed, a sizable relic abundance of weakly interacting massive par- 
ticles (WIMPs) is generally expected to be produced as a byproduct of 
our universe's hot youth. In this section, I discuss this process and the 
determination of the relic abundance of a WIMP [13, 14, 15]. 

Consider a stable particle, X, which interacts with Standard Model 
particles, Y, through some process XX YY (or XX YY if X is its 
own antiparticle). In the very early universe, when the temperature was 
much higher than nix , the processes of XX creation and annihilation were 
equally efficient, leading X to be present in large quantities alongside the 
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various particle species of the Standard Model. As the temperature of the 

universe dropped below mx, however, the process of XX creation became 
exponentially suppressed, while XX annihilation continued unabated. In 
thermal equilibrium, the number density of such particles is given by 



where gx is the number of internal degrees of freedom of X. 

If these particles were to remain in thermal equilibrium indefinitely, 

their number density would become increasingly suppressed as the universe 
cooled, quickly becoming cosmo logically irrelevant. There are ways that 
a particle species might hope to avoid this fate, however. For example, 
baryons are present in the universe today because of a small asymmetry 
which initially existed bctwcx^n tlic^ number of baryons and antibaryons; 
when all of the antibaryons had annihilated with baryons, a small residual 
of baryons remained. The baryon-antibaryon asymmetry prevented the 
complete annihilation of these particles from taking place. 

While it is possible that a particle-antiparticle asymmetry is also behind 
the existence of dark matter, there is an even simpler mechanism which can 
lead to the survival of a sizable relic density of weakly interacting particles. 
In particular, the self-annihilation of weakly interacting species can be con- 
tained by the competing effect of Hubble expansion. As the expansion and 
corresponding dilution of WIMPs increasingly dominates over the annihi- 
lation rate, the number density of X particles becomes sufficiently small 
that they cease to interact with each other, and thus survive to the present 
day. Quantitatively, the competing effects of expansion and annihilation 
are described by the Boltzmann equation: 



where nx is the number density of WIMPs, H = R/R = (Stt^p/SMpO^/^ jg 
the expansion rate of the universe, and < ajf^jf | > is the thermally aver- 
aged XX annihilation cross section (multiplied by their relative velocity). 

Prom Eq. 5, we can identify two clear limits. As I said before, at very 
high temperatures (T ^ mx) the density of WIMPs is given by the equi- 
librium value, nx, eq- In the opposite limit {T <C mx), the equilibrium 
density is very small, leaving the terms 3Hnx and < •J^xl'^l > ''^x each 
further deplete the number density. For sufficiently small values of nx , the 
annihilation term becomes insignificant compared to the dilution due to 





dnx 
dt 



+ 3Hnx = - < (7xx\v\ > {nx - '^X,eq)> 



(5) 
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Figure 1. A schematic of the comoving number density of a stable species as it evolves 
through the process of thermal freeze-out. 



Hubble expansion. When this takes place, the comoving number density of 
WIMPs becomes fixed — thermal freeze-out has occurred. 

The temperature at which the number density of the species X departs 
from equilibrium and freezes out is found by numerically solving the Boltz- 
mann equation. Introducing the variable x = mx/T, the temperature at 
which freeze-out occurs is approximately given by 



mx , 
a;FO = 7^ « In 
-tpo 



[45 gx mxMpi{a + 6b/xFo) 

''^''■^ 8 27r3 1/21/2 
9* sjpo 



(6) 



Here, c ^ 0.5 is a numerically determined quantity, gf* is the number of 

external degrees of freedom available (in the Standard Model, g^, ^ 120 at 
T ~ 1 TeV and (7* ~ 65 at T ~ 1 GeV) , and a and b are terms in the non- 
relativistic expansion, < ctj^xI'*^! >= a + b < v"^ > +0{v^). The resulting 
density of WIMPs remaining in the universe today is approximately given 

by 

2 1.04 X 109 GeV-^ xpo 

iixii ~ — -rTT, . in 

Mm gll\a + mix^o) 



5 



If X has a GeV-TeV scale mass and a roughly weak-scale annihilation cross 
section, freeze-out occurs at x-po ~ 20 — 30, resulting in a relic abundance 
of 

In other words, if a GeV-TeV scale particle is to be thermally produced with 
an abundance similar to the measured density of dark matter, it must have 
a thermally averaged annihilation cross section on the order of 3 x 10^^^ 
cm^/s. Remarkably, this is very similar to the numerical value arrived 
at for a generic weak-scale interaction. In particular, a^/(100GeV)^ ~ 
pb, which in our choice of units (and including a factor of velocity) is 
3 X 10~^^ cm^/s. The similarly between this result and the value required 
to generate the observed quantity of dark matter has been dubbed the 
"WIMP miracle". While far from constituting a proof, this argument has 
lead many to conclude that dark matter is likely to consist of particles 
with weak-scale masses and interactions and certainly provides us with 
motivation for exploring an electroweak origin of our universe's missing 
matter. 

2.1. Case Example - The Thermal Abundance of a Light or 
Heavy neutrino 

At first glance, the Standard Model itself appears to contain a plausible 
candidate for dark matter in the form of neutrinos. Being stable and weakly 
interacting, neutrinos are a natural place to start in our hunt for dark 
matter's identity. 

In the case of a Standard Model neutrino species, the relatively small 
annihilation cross section (< a\v\ >~ lO""^^ cm^/s) and light mass leads to 
an overall freeze-out temperature on the order of Tpo ^ MeV, and a relic 
density of 

O.+./.^«0.l(^). (9) 

As constraints on Standard Model neutrino masses require to be well 
below 9 eV, we are forced to conclude that only a small fraction of the dark 
matter could possibly consist of Standard Model neutrinos. Furthermore, 
even if these constraints did not exist, such light neutrinos would be highly 
relativistic at the time of freeze-out {T^o/mv ^ MeV/m,^ ^ 1) and thus 
would constitute hot dark matter, in conflict with observations of large 
scale structure. 
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Moving beyond the Standard Model, we could instead consider a heavy 

4th generation Dirac neutrino. In this case, the annihilation cross section 
can be much larger, growing with the square of the neutrino's mass up 
to the Z pole, rriv ~ mz/2, and declining with above ^ mz/2. 
For a GcV-TcV mass neutrino, the process of freeze-out yields a cold relic 
{T-po/^v ~ 0(0. 1)), with an abundance approximately given by 

—] , MeV<m^<mz/2 (10) 

nil, J 

n^+.h? ~ 0.1 {j^^^ , « m.. (11) 

Thus, if the relic density of a heavy neutrino species is to constitute the 
bulk of the observed dark matter abundance, we find that it must have a 

mass of approximately 5 GeV [16], or several hundred GcV. The former 
case is excluded by LEP's measurement of the invisible width of the Z, 
however, which rules out the possibility of a fourth active neutrino species 
lighter than half of the Z mass [17]. The later case, although consistent 
with the bounds of LEP and other accelerator experiments, is excluded by 
the limits placed by direct dark matter detection experiments (which I will 
discuss in Sec. 4). 



2.2. A Detour Into Coannihilations 

In some cases, particles other than the WIMP itself can play an important 
role in the freeze-out process [18]. Before such a particle can significantly 
impact the relic density of a WIMP, however, it must first manage to be 
present at the temperature of freeze-out. 

The relative abundances of two species at freeze-out can be very roughly 
estimated by 

^ ^ ^""^f ^° . (12) 

Considering, for example, a particle with a mass twice that of the WIMP 
and a typical freeze-out temperature of mx/T^o ~ 20, there will be only 
~ e~'"^/e^^° ^ 10^^ Y particles for every X at freeze-out, thus making Y 
completely irrelevant. If my were only 10% larger than mx, however, we 
estimate ny jnx ~ e~^^/e~^° 10~^. In this quasi-degenerate case, the 
additional particle species can potentially have a significant impact on the 
dark matter relic abundance. 
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To quantitatively account for other species in the calculation of the relic 
abundance of a WIMP, we make the following substitution (for both a and 
b) into Eqs. 6 and 7: 

<7An„ - aM^) =I^^i,j^^i^ + + A,)3/2e-(^.+^.), (13) 

where the double sum is over all particle species {i,j=l denoting the WIMP 
itself) and aij is the cross section for the coannihilation of species i and j (or 
self-annihilation in the case of z = j) into Standard Model particles. As the 
eifective annihilation cross section has a strong dependence on a;, we must 
integrate Eqs. 6 and 7 over x (or T). The quantities Aj = (m^ — mi) /mi 
denote the fractional mass splittings between the species i and the WIMP. 
The effective number of degrees of freedom, gEs{x), is given by: 

gM^) = Y.9i{^ + Ai)'/'e--^'. (14) 

i 

To better understand how the introduction of particles other than the 
WIMP can effect the process of freeze-out, lets consider a few simple cases. 
First, consider one additional particle with a mass only slightly above that 
of the WIMPs (A2 <^ 1), and with a comparatively large coannihilation 
cross section, such that ai^2 ^ this case, gES ^ gi + 52, and 

(TEff « 0-1,2.91.92/(51 + <72)^- Since (t-e« is much larger than the WIMP's 
self-annihilation cross section, the relic density of WIMPs will be sharply 
suppressed. This is the case that is usually meant by the term "coannihi- 
lation" . 

Alternatively, consider the opposite case in which the WIMP and the 
additional quasi-degenerate particle do not coannihilate efficiently (ui 2 «C 
0-1,1, 0-2,2)- Here, (Jes « 0-1,1 .g^/(.gi+,g2)^ + 02, 2 52/(51+52)^, which in some 
cases can actually be smaller than that for the process of self-annihilation 
alone, leading to an enhanced relic abundance. Physically speaking, what 
is going on here is that the two species are each freezing out independently 
of each other, after which the heavier species decays, producing additional 
WIMPs as a byproduct. 

As an extreme version of this second case, consider a scenario in which 
the lightest state is not a WIMP, but is instead a purely gravitationally 
interacting particle. A slightly heavier particle with weak interactions will 
self-annihilate much more efficiently than it will coannihilate with the light- 
est particle (0-1,2 is negligible), leading the two states to freeze-out inde- 
pendently. The gravitationally interacting particle, however, never reaches 
thermal equilibrium, so could potentially have not been produced in any 
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significant quantities up until this point. Well after freezing out, the heavier 

particles will eventually decay, producing the stable gravitationally inter- 
acting lightest state. Although the resulting particles are not WIMPs (they 
do not have weak interactions), they are naturally produced with approx- 
imately the measured dark matter abundance because of the WIMP-like 
properties of the heavier state. In other words, this case - known as the 
"superWIMP" scenario [19] - makes use of the coincidence between the 
electroweak scale and the measured dark matter abundance without the 
dark matter actually consisting of WIMPs. Because gravitationally inter- 
acting particles and other much less than weakly interacting particles are 
almost impossible to detect astrophysically, superWIMPs are among the 
dark matter hunter's worst nightmares. 

3. Beyond The Standard Model Candidates For Dark 
Matter 

There has been no shortage of dark matter candidates proposed over the 
years. A huge variety of beyond the Standard Model physics models have 
been constructed which include a stable, electrically neutral, and colorless 
particle, many of which could serve as a phenomenologically viable candi- 
date for dark matter. I could not possibly list, must less review, all of the 
proposed candidates here. Finding the "WIMP miracle" (as discussed in 
Sec. 2) to be fairly compelling (along with the hierarchy problem, which 
strongly suggests the existence of new particles at or around the electroweak 
scale), I chose to focus my attention on dark matter in the form of weak- 
scale particles. So although the dark matter of our universe could plausibly 
consist of particles ranging from 10^^ eV axions to 10^^ GeV WIMPzillas, I 
will ignore everything but those particle physics frameworks which predict 
the existence of a stable, weakly interacting particle with a mass in the few 
GeV to few TeV range. 

3.1. Supersymmetry 

For a number of reasons, supersymmetry is considered by many to be among 
the most attractive extensions of the Standard Model. In particular, weak- 
scale supersymmetry provides us with an elegant solution to the hierarchy 
problem [20] , and enables grand unification by causing the gauge couplings 
of the Standard Model to evolve to a common scale [21]. From the stand- 
point of dark matter, the lightest superpartner is naturally stable in models 
that conserve i?-parity. i?-parity is defined as i? = (_ 1)35-1-^+25 (-^^ g^j^^ 
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S denoting baryon number, lepton number and spin), and thus is assigned 
as i? = +1 for all Standard Model particles and 7? = —1 for all supcrpart- 
ners. i?-parity conservation, therefore, requires superpartners to be created 
or destroyed in pairs, leading the lightest supersymmetric particle (LSP) to 
be stable, even over cosmological timescales. 

The identity of the LSP depends on the hierarchy of the STipersymmetric 
spectrum, which in turn is determined by the details of how supersymmetry 
is broken. The list of potential LSPs which could constitute a plausible dark 
matter candidate is somewhat short, however. The only electrically neutral 
and colorless superparnters in the minimal supersymmetric standard model 
(MSSM) are the four neutralinos (superpartners of the neutral gauge and 
Higgs bosons), three sneutrinos, and the gravitino. The lightest neutralino, 
in particular, is a very attractive and throughly studied candidate for dark 
matter [22]. 

Before discussing the phenomenology of neutralino dark matter, lets 
briefly contemplate the possibility that sneutrinos might make up the; dark 
matter of our universe. In many respects, sneutrino dark matter would 
behave very similarly to a heavy 4th generation neutrino, as discussed in 
Sec. 2.1. In particular, like neutrinos, sneutrinos are predicted to annihilate 
to Standard Model fermions efficiently through the s-channel exchange of 
a Z boson (as well as through other diagrams). As a result, sneutrinos 
lighter than about 500-1000 GeV would be under produced in the early 
universe (a ~10 GeV sneutrino would also be produced with approximately 
the measured dark matter abundance, but is ruled out by LEP's invisible 
Z measurement). 

The Feynman diagram corresponding to sneutrino annihilation into 
quarks through a s-channel Z exchange can be turned on its side to pro- 
duce an elastic scattering diagram with quarks in nuclei. When the elastic 
scattering cross section of a ~100-1000 GeV sneutrino is calculated, we 
find that it is several orders of magnitude larger than current experimen- 
tal constraints [23]. We are thus forced to abandon MSSM sneutrinos as 
candidates for dark matter. 

In the MSSM, the superpartners of the four Standard Model neutral 
bosons (the bino, wino and two neutral higgsinos) mix into four physical 
states known as neutralinos. Often times, the lightest of these four states 
is simply referred to as "the neutralino" . The neutralino mass matrix can 
be used to determine the masses and mixings of these four states. In the 
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B-W^-Hi-H2 basis, this matrix is given by 



= 



(15) 




-mz cos P sin 9 sin /? sin \ 
TOz cos /3 cos —mzsmPcosOyy 



—mz cos I3sm6^r mz cos P cos 0^ 
\ mz sin /? sin —mz sin (3 cos 6^ 



'w 



'w 




where Mi and M2 are the bino and wino masses, fj, is the higgsino mass 
parameter, 6w is the Weinberg angle, and tan/3 = vijvi is the ratio of 
the vacuum expectation values of the Higgs doublets. This matrix can be 
diagonalized into mass eigenstates by the unitary matrix A'', 



In terms of the elements of the matrix, A/', the lightest neutralino is given 
by the following mixture of gaugino and higgsino components: 



The quantities [iVnp + |iVi2p and |A^i3p + [iV^p are often referred to 
as the gaugino fraction and higgsino fraction of the lightest neutralino, 
respectively. 

The lightest neutralino can annihilate through a wide variety of Feyn- 
man diagrams. In Fig. 2, we show some of the most important of these; 
although it is far from an inclusive list. Which of these diagrams dominate 
the process of thermal freeze-out in the early universe depends on the com- 
position of the lightest neutralino, and on the masses and mixings of the 
exchanged particles [24] . Since so many different diagrams can potentially 
contribute to neutralino annihilation (not to mention the many possible 
coannihilation processes [25, 26, 27, 28]), the resulting relic density de- 
pends on a large number of supersymmetric parameters and is not trivial 
to calculate accurately. Publicly available tools such as DarkSUSY [29] and 
MicroOmegas [30] are often used for this purpose. 

The mass and composition of the lightest neutralino is a function of 
four supersymmetric parameters: Mi, M2, \x and tan/3. This becomes 
further simplified if the gaugino masses are assumed to evolve to a single 
value at the GUT scale, yielding a ratio at the electroweak scale of Mi = 
I tan^ 6wM2 ~ O.5M2. In this case, the lightest neutralino has only a 
small wino fraction and is largely bino-like (higgsino-like) for Mi <C \fi\ 



(16) 



+ N13H1 + N14H2. 



(17) 



(Ml » ImI). 
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Over much or most of the supersymmetric parameter space, the relic 

abundance of neutralinos is predicted to be in excess of the observed dark 
matter density. To avoid this, we are forced to consider the regions of 
parameter space which lead to especially efficient neutralino annihilation in 
the early universe. In particular, the following scenarios are among those 
which can lead to a phenomenologically viable density of neutralino dark 
matter: 



• If the lightest neutralino has a significant higgsino or wino fraction, 
it can have fairly large couplings and, as a result, annihilate very 
efficiently. 

• If the mass of the lightest neutralino is near a resonance, such as the 
CP-odd Higgs pole, it can annihilate efficiently, even with relatively 
small couplings. 

• If the lightest neutralino is only slightly lighter than another su- 
perpartner, such as the lightest stau, coannihilations between these 
two states can very efficiently deplete the dark matter abundance. 



In Fig. 3, we illustrate these regions within the context of a specific 
subset of the MSSM known as the CMSSM (C stands for constrained). In 
this framework, all of the scalar masses are set to a common value mo at the 
GUT scale, from which the electrowcak scale values are determined by RGE 
evolution. Similarly the three gaugino masses are each set to mi/2 at the 
GUT scale. In each frame, the narrow blue regions denote the parameter 
space in which neutralino dark matter is predicted to be generated with the 
desired abundance (0.0913 < fi^o/i^ < 0.1285). In the corridor along side 
of the LEP chargino bound (m^± > 104 GeV), ^ and Mi are comparable 
in magnitude, leading to a mixed bino-higgsino LSP with large couplings. 
Within the context of the CMSSM, this is often called the "focus point" 
region. In the bottom portion of each frame, the lightest stau (f) is the LSP, 
and thus does not provide a viable dark matter candidate. Just outside of 
this region, however, the stau is slightly heavier than the lightest neutralino, 
leading to a neutralino LSP which efficiently coannihilates with the nearly 
degenerate stau. In the lower right frame, a viable region also appears 
along the CP-odd Higgs resonance {m-^o w niA/^)- This is often called the 
A-funnel region. 
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Figure 3. Representative regions of the CMSSM parameter space. The blue regions 
predict a neutrahno density consistent with the measured dark matter abundance. The 
shaded regions to the upper left and lower right are disfavored by the LEP chargino 
bound and as a result of containing a stau LSP, respectively. The LEP bound on the 
light Higgs mass is shown as a solid line (m/j = 114 GeV). The region favored by 
measurements of the muon's magnetic moment are shown as a light shaded region (at 
the 3a confidence level) [31]. In each frame, we have used Aq = and fi > 0. (Figures 
generated by Gabriel Caceres, using the package DarkSUSY [29].) 

3.2. Kaluza-Klein Dark Matter in Models With Universal 
Extra Dimensions 

Supersymmetry is not by any means the only particle physics framework 
from which viable a dark matter candidate can arise. As an alternative, I 
will discuss in this section the possibility of Kaluza-Klein dark matter in 
models with a universal extra dimension. 

In recent years, interest in theories with extra spatial dimensions has 
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surged. In particular, a colossal amount of attention has been given to 

two classes of extra dimensional theories over the past decade: scenar- 
ios featuring one or more large (millimeter-scale), flat extra dimensions 
(ADD) [32, 33], and the Randall- Sundrum scenario, which introduces an 
additional small dimension with a large degree of spatial curvature [34] . A 
somewhat less studied but interesting class of extra-dimensional models, 
which goes by the name of universal extra dimensions (UED), postulates 
the existence of a flat extra dimension (or dimensions) through which all 
of the Standard Model (SM) fields arc; free to propagate (rather than being 
confined to a brane as some or all of them are in the ADD and Randall- 
Sundrum models) [35]. 

For a number of reasons, extra dimensions of size R ~ TeV~^ are par- 
ticularly well motivated within the context of UED [36] (much smaller than 
those found in the ADD model). Among these reasons is the fact that a 
TeV-scale Kaluza-Klein (KK) state, if stable, colorless, and electrically neu- 
tral, could potentially serve as a viable candidate for dark matter [37, 38] . 

Standard Model fields with momentum in an extra dimension appear 
as heavy particles, called KK states. This leads to a tower of KK states for 
each Standard Model field, with tree-level masses given by: 

™xw = ^ +"^x(o). (18) 

where X^") is the nth Kaluza-Klein excitation of the Standard Model field, 
X, and i? ~ TeV"^ is the size of the extra dimension. X^"^ denotes the 
ordinary Standard Model particle (known as the zero mode). 

If the extra dimensions were simply wrapped (compactified) around a 
circle or torus, then extra dimensional momentum conservation would en- 
sure the conservation of KK number (n) and make the lightest first level 
KK state stable. Realistic models, however, require an orbifold to be intro- 
duced, which leads to the violation of KK number conservation. A remnant 
of KK number conservation called KK-parity. however, can remain and lead 
to the stability of the lightest KK particle (LKP) in much the same way 
that J?-parity conservation prohibits the decay of the lightest supersym- 
metric particle. 

In order for the LKP to be a viable dark matter candidate, it must be 
electrically neutral and colorless. Possibilities for such a state include the 
first KK excitation of the photon, Z, neutrinos, Higgs boson, or graviton. 
Assuming that is considerably larger than any of the Standard Model 
zero mode masses, Eq. 18 leads us to expect a highly degenerate spectrum 



15 



of Kaluza-Klein states at each level (although this picture is somewhat 
modified when radiative corrections and boundary terms are included). Of 
our possible choices for the LKP, the relatively large zero-mode mass of the 
Higgs make its first level KK excitation an unlikely candidate. Further- 
more, KK neutrinos are excluded by direct detection experiments, just as 
sneutrinos or heavy 4th generation Dirac neutrinos are. For these reasons, 
we focus on the mixtures of the KK photon and KK Z as our dark matter 
candidate (note that, unlike higgsino and gauginos, the KK Higgs has a 
different spin than the KK photon and KK Z, and thus does not mix with 
these states). 

The mass eigenstates of the KK photon and KK Z are very nearly 
identical to their gauge eigenstates, B^"^ and The reason for this 

can be seen from their mass matrix: 

/ ^ + + \g\v'' ^ \g^g2V^ \ ^^^^ 

Here « 174 GeV is the Higgs vacuum expectation value. In the well 
known zero- mode case (n = 0), there is significant mixing between B^^^ and 
p[/3(o) ^gjj^2 0^ ^ 0.23). In the absence of radiative corrections (5m^(„) = 
(5to^(„) = 0), the same mixing angle is found at the first KK level as well. If 
the difference between (5m^(„) and 5to^(„) is larger than the (rather small) 
off diagonal terms, however, the mixing angle between these two KK states 
is driven toward zero. Using typical estimates of these radiative corrections, 
the effective first KK level Weinberg angle is found to be approximately 
siia^ 0w,i ~ 10~^. Thus the mass eigenstate often called the "KK photon" 
is not particularly photon-like, but instead is nearly identical to the state 

The KK state S'^' annihilates largely to Standard Model (zero-mode) 
fermions through the t-channel exchange of KK fermions, with a cross sec- 
tion given by 

^ 32,256^ 7rm^(i, ^ ' 

As the B^^^ couples to the fermions' hypercharge, the cross section scales 
with -|- and most of its annihilations proceed to charged lepton 
pairs. 

If Fig. 4, the thermal relic abundance of KK dark matter is shown as 
a function of its mass. Because the first level KK spectrum is expected to 
be quasi-degenerate, coannihilations are likely to play an important role. 
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Figure 4. The thermal reUc abundance of KK dark matter (f?'^') without the effects 
of any other KK species (solid line) and including the effects of KK leptons 5% and 1% 
heavier than the LKP (dashed and dotted linos). Shown as a horizontal band is the 
measured dark matter abundance [6]. Adapted from Ref. [37]. 

In the figure, results are shown ignoring the effects of other KK states 
(solid line) and including the effects of KK leptons 5% or 1% heavier than 
the LKP (dashed and dotted lines, respectively). Note that the KK lep- 
tons lead to a larger relic abundance, due to the fact that they freeze-out 
quasi-independently from the LKP and then increase the number of LKPs 
through their decays. Depending on the details of the KK spectrum, many 
other states could potentially effect the relic density of KK dark matter 
as well. LKP masses from approximately 500 GeV to several TeV can po- 
tentially lead to a relic density consistent with the measured dark matter 
abundance [39, 40]. 

3.3. A Note On Other Possibilities For TeV-Scale Dark 
Matter 

At this time, I would like to make a general comment about some of the 
many other possibilities for the particle identity of dark matter. It is in- 
teresting to note that a wide range of solutions to the gauge hierarchy 
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problem also introduce a candidate for dark matter. In particular, in order 

to be consistent with electrowcak precision measurements, new physics at 
the TeV scale typically must possess a discrete symmetry. Supersymmetry, 
for example, accomplishes this through i?,-parity conservation. Similarly, 
Little Higgs models can avoid problems with electrowcak precision tests by 
introducing a discrete symmetry called T-parity. And just as _R-parity sta- 
bilizes the lightest superpartner, T-parity can enable the possibility of dark 
matter by stabilizing the lightest T-parity odd state in the theory [41]. So, 
thinking outside of any specific particle physics framework, we can imag- 
ine a much larger class of TeV-scale physics scenarios which address the 
hierarchy problem without violating electrowcak precision measurements 
by introducing a discrete symmetry, which in turn leads to a stable dark 
matter candidate. In this way, many of the motivations and much of the 
phenomenology discussed within the context of supersymmetry can actually 
be applied to a more general collection of particle physics models associated 
with the electrowcak scale. 



4. Direct Detection 

Turning our attention now to dark matter detection, we begin with those ex- 
periments which attempt to detect dark matter particles through their elas- 
tic scattering with nuclei, including CDMS [42], XENON [43], ZEPLIN [44], 
EDELWEISS [45], CRESST [46], CoGeNT [47], DAM A/LIBRA [48], 
COUPP [49], WARP [50], and KIMS [51]. This class of techniques is col- 
lectively known as direct detection, in contrast to indirect detection efforts 
which attempt to observe the annihilation products of dark matter parti- 
cles. 

A WIMP striking a nucleus will induce a recoil of energy given by 

^ ^ ^ |gT ^ 2/X^f^(l-C0Sg) ^ m^^MnuclensvH^-COsO) ^21) 
2Mnucleus 2Mnucleus {mx + Afnucleus)^ 

where q is the WIMP's momentum, v is its velocity, and fj, is the reduced 
mass. For mx ^ ^nucleus and a velocity of ^300 km/s, we expect typical 
recoil energies of E'recoii ~ -Mnucieust"^ ~ 1-100 keV. 

WIMPs scatter with nuclei in a target at a rate given by 

JBmin Jvrnir, TUX d\q\ ^ ' 

where p is the dark matter density, a is the WIMP-nuclei clastic scattering 
cross section, and f{v) is the velocity distribution of WIMPs. The limits 



18 



of integration are set by the galactic escape velocity, t;max ~ 650 km/s, and 
kincmatically by Wmin = (-Erocoii The minimum energy is 
set by the energy threshold of the detector, which is typically in the range 
of several keV to several tens of keV. 

WIMPs can potentially scatter with nuclei through both spin- 
independent and spin-dependent interactions. The experimental sensitiv- 
ity to spin-independent couplings benefits from coherent scattering, which 
leads to cross sections (and rates) proportional to the square of the atomic 
mass of the target nuclei. The cross sections for spin-dependent scattering, 
in contrast, are proportional to J(J -|- 1), where J is the spin of the target 
nucleus, and thus do not benefit from large target nuclei. As a result, the 
current experimental sensitivity to spin-dependent scattering is far below 
that of spin-independent interactions. For this reason, wc consider first the 
case of spin-independent scattering of WIMPs with nuclei (we will return 
to spin-dependent scattering in Sec. 5.3). 

The spin-independent WIMP-nucleus elastic scattering cross section is 
given by 

4?77 ^ 777^ 
Cr « X° nucleus ^ 

7r(m^o + mnucleus) 

where Z and A are the atomic number and atomic mass of the nucleus, fp 
and /„ are the WIMP's couplings to protons and neutrons, given by [24] 

q=u,d^s q=c,b,t 

where a, are the WIMP-quark couplings and /^^^ » 0.020 ± 0.004, f^^J w 
0.026 ± 0.005, f^^ w 0.118 ± 0.062, f^^ w 0.014 ± 0.003, f^^ w 0.036 ± 
0.008, /^"^ « 0.118 ± 0.062 are quantities measured in nuclear physics 
experiments [52]. The first term in Eq. 24 corresponds to interactions with 
the quarks in the target nuclei. The second term corresponds to interactions 
with the gluons in the target through a colored loop diag ram. is given 
by 1 - - f^^J - f^^^ « 0.84, and analogously, f^"^ « 0.83. 

4.1. Direct Detection of Neutralino Dark Matter 

Neutralinos can elastically scatter with quarks through either f-channel 
CP-even Higgs exchange, or s-channel squark exchange: 
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In addition to these diagrams, we can write analogous processes in which 
the WIMP couples to giuons in the target through a quark/squark loop. By 
calculating the WIMP-quark couplings, a^, we can also implicitly include 
the interactions of neutralinos with giuons in the target nuclei as well (see 
Eq. 24). 

The neutralino-quark coupling, in which all of the supersymmetry 
model-dependent information is contained, is given by [53] 
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where throughout z = 1 for up-type quarks and z = 2 for down type quarks. 
fn\i,m2i denote the squark mass eigenvalues and rj is the matrix which 
diagonalizes the squark mass matrices, diag{m\,m2) = r]M^ri~^. yi, Tsi 
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and Cj denote hypercharge, isospin and electric charge of the quarks. For 
scattering off of up-type quarks 

Si = Nis, 82 = Nu, B = smP, C = sina, i:» = cosQ!, (27) 

whereas for down-type quarks 

61 = Ni4, 52 = —Ni3, B = cos/3, C = cosQ!, £> = — sina. (28) 

The quantity a is the angle that diagonalizes the CP-even Higgs mass 
matrix. 

The first two terms in Eq. 25 correspond to interactions through the 
exchange of a squark, while the final term is generated through Higgs ex- 
change. To help develop some intuition for what size neutralino-nucleon 
cross sections we might expect, lets consider a few simple limits. First, 
consider the case in which the scattering is dominated by heavy Higgs (H) 
exchange through its couplings to strange and bottom quarks. This behav- 
ior is often found when the squarks are heavy and cos q w 1 (which implies 
moderate to large tan/3 and niA ^ mH ~ TOij±)- In this case, the leading 
contribution to the neutralino-nucleon cross section is 



possibility of very large rates. For a ^100 GeV neutralino, for example, 
a 200 GeV heavy Higgs mass leads to cross sections with nucleons on the 
order of IQ-^ to 10"^ pb for ~ 200 GeV, or 10"^ to IQ-^ pb for \n\ ~ 1 
TeV. 

Second, we can consider the case in which the cross section is dominated 
by light Higgs boson (h) exchange through its couplings to up-type quarks. 
This is often found in the case of heavy squarks and heavy to moderate H. 
In this limiting case 



If the heavy Higgs {H) is heavier than about -^500 GeV, exchange of the 
light Higgs generally dominates, leading to cross sections of around 10~^ 
to 10-1° pb for in the range of 200 GeV to 1 TeV. 

Third, consider the case in which the elastic scattering cross section is 
dominated by the exchange of squarks through their couplings to strange 
and bottom quarks. This is found for large to moderate tan/3 and squarks 




Here, the cross section scales with rrif/ and with tan^/3, a leads to the 
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with masses well below 1 TeV. In this limiting case, and with approximately 
diagonal squark mass matrices, 

"^"^ 4.m^cos^/3m| i,^"- + 27^ ' 'iommated,tan/3 » 1). 

(31) 

For squarks lighter than ~1 TeV, squark exchange can potentially provide 
the dominant contribution to neutralino-nuclei elastic scattering. 

Using Eq. 22, we can crudely estimate the minimum target mass re- 
quired to potentially detect neutralino dark matter. A detector made up 
of Germanium targets (such as CDMS or Edelweiss, for example) would 
expect a WIMP with a nucleon-level cross section of 10~^ pb (10~^^ cm^) 
to yield approximately 1 elastic scattering event per kilogram-day of ex- 
posure. Such a target mass could thus be potentially sensitive to strongly 
mixed gaugino-higgsino neutralinos with light mjj and large tan/?. The 
strongest current limits on spin-independent scattering have been obtained 
using ~10^ kilogram-days of exposure. In contrast, reaching sensitivities 
near the 10^^° pb level will require ton-scale detectors capable of operating 
for weeks, months or longer with very low backgrounds. 

In Fig. 5, wc show the current constraints on the WIMP-nucleon spin- 
independent elastic scattering cross section, as a function of the WIMP's 
mass. The most stringent constraints currently come from the CDMS and 
XENON-10 collaborations, which each have obtained limits at around the 
10"'^ pb level. The CDMS and XENON-100 collaborations are each an- 
ticipated to place limits a factor of several times more stringent within 
the next year (unless a positive detection is made, that is). Although the 
more distant future is more difficult to project, it is generally expected that 
experiments approaching the ton-scale will reach sensitivies near 10~^ or 
10~^° pb within in the next few to several years. 

4.2. Direct Detection of Kaluza-Klein Dark Matter 

In the case of Kaluza-Klein dark matter (as described in Sec. 3.2), the 
WIMP-quark coupling a, receives contributions from the s-channel ex- 
change of KK quarks and the t-channel Higgs boson exchange [38, 54]. 
This leads to a contribution to the cross section from Higgs exchange 
which is proportional to l/(m^(i) m^) and a contribution from KK-quark 
exchange which is approximately proportional to l/(m^(i, A'*), where 
A = (mq(i) — ms(i))/mB(i) is the fractional mass splitting of the KK quarks 
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WIMP Mass [GeV/c^] 

Figure 5. Current constraints on the WIMP-nucleon, spin-independent elastic scatter- 
ing cross section. From bottom-to-top on the right side of the figure, the lines corre- 
spond to the limits from the CDMS [42], XENON-10 [43], WARP [50], CRESST [46], 
ZEPLIN [44], and Edelweiss [45] experiments. This figure was generated using the Dark 
Matter Limit Plotter [55]. 
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With such small cross sections, we will most likely have to wait for 
ton-scale detectors before this model will be tested by direct detection ex- 
periments. 
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4.3. Some Model Independent Comments Regarding Direct 
Detection 

In the cases of the two dark matter candidates discussed thus far in this 
section, we are lead to expect a rather smah elastic scattering cross sec- 
tions between WIMPs and nuclei typically below or well below current 
experimental constraints. This is not a universal prediction for a generic 
WIMP, however. In fact, both neutralinos and Kaluza-Klein dark matter 
represent somewhat special cases in which direct detection rates are found 
to be particularly low. 

To illustrate this point, consider a Dirac fermion or a scalar WIMP 
which annihilates in the early universe to fermions with roughly equal cou- 
plings to each species - a heavy 4th generation neutrino or sneutrino, for 
example. We can take the Feynman diagram for the process of this WIMP 
annihilating to quarks and turn it on its side, and then calculate the result- 
ing elastic scattering cross section. What we find is that, if the interaction 
is of scalar or vector form, such a WIMP will scatter with nuclei several 
orders of magnitude more often than is allowed by the limits of CDMS, 
XENON and other direct detection experiments. Similar conclusions are 
reached for many otherwise acceptable WIMP candidates [56]. A warning 
well worth keeping in mind for any WIMP model builder is, "Beware the 
crossing symmetry!" . 

So what is it about neutralinos or Kaluza-Klein dark matter than enable 
them to evade these constraints? In the case of neutralinos, the single most 
important feature is the suppression of its couplings to light fermions. Being 
a Majorana fermion, a neutralino's annihilation cross section to fermion 
pairs (at low velocity) scales with av cx TvPj/m^o. As a result, neutralinos 
annihilate preferentially to heavy fermions (top quarks, bottom quarks, 
and taus) or gauge/Higgs bosons. As heavy fermions (and gauge/Higgs 
bosons) are largely absent from nuclei, the potentially dangerous crossing 
symmetry does not apply. More generally speaking, current direct detection 
constraints can be fairly easily evaded for any WIMP which interacts with 
quarks through Higgs exchange, as the Yukawa couplings scale with the 
fermion's mass. 

Alternatively, if the WIMP's couplings are simply very small, direct 
detection constraints can also be evaded. Small couplings, however, leave 
us in need of a mechanism for efficiently depleting the WIMP in the early 
universe. But even with very small couplings, a WIMP might efficiently 
coannihilate in the early universe, or annihilate through a resonance, lead- 
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ing to an acceptable relic abundance. In this way, coannihilations and 
resonances can considerably suppress the rates expected in direct detection 
experiments. 

5. Indirect Detection 

Direct detection experiments are not the only technique being pursued in 
the hope of identifying the particle nature of dark matter. Another major 
class of dark matter searches are those which attempt to detect the products 
of WIMP annihilations, including gamma rays, neutrinos, positrons, elec- 
trons, and antiprotons. These methods are collectively known as indirect 
detection and are the topic of this section. 

5.1. Gamma Rays From WIMP Annihilations 

Searches for the photons generated in dark matter annihilations have a 
key advantage over other indirect detection techniques in that these par- 
ticles travel essentially unimpeded. In particular, unlike charged particles, 
gamma rays are not deflected by magnetic fields, and thus can potentially 
provide valuable angular information. For example, point-like sources of 
dark matter annihilation radiation might appear from high density regions 
such as the Galactic Center or dwarf spheroidal galaxies. Furthermore, over 
galactic distance scales, gamma rays are not attenuated, and thus retain 
their spectral information. In other words, the spectrum that is measured 
is the same as the spectrum generated in the dark matter annihilations. 

The spectrum of photons produced in dark matter annihilations de- 
pends on the details of the WIMP being considered. Neutralinos, for ex- 
ample, typically annihilate to final states consisting of heavy fermions (bb, 
tt, T+T-) or gauge and/or Higgs bosons (ZZ, W+W^ , HA, hA, ZH, Zh, 
ZA, W^H^, where H, h, A and are the Higgs bosons of the minimal 
supersymmetric standard model) [24]. With the exception of the r+r" 
channel, each of these annihilation modes result in a very similar spectrum 
of gamma rays. The gamma ray spectrum from a WIMP which annihilates 
to leptons can be quite different, however. This can be particularly impor- 
tant in the case of Kaluza-Klein dark matter in models with one universal 
extra dimension, for example, in which dark matter particles annihilate 
significantly to e+e", /x^/U" and T+r~ [37, 38]. In Fig. 6, we show the pre- 
dicted gamma ray spectrum, per annihilation, for several possible WIMP 
annihilation modes. 

In addition to generating continuum gamma rays through the decays of 
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Figure 6. The gamma ray spectrum per WIMP annihilation for a 100 GeV (left) and 500 
GeV (right) WIMP. Each curve denotes a different choice of the dominant annihilation 
mode: bb (solid cyan), ZZ (magenta dot-dashed), W'^W~ (blue dashed), t+t~ (black 
solid), e~^e~ (green dotted) and fJ.^ li~ (red dashed). 



quarks, leptons, Higgs bosons or gauge bosons, dark matter particles can 
produce gamma rays directly, leading to monoenergetic spectral signatures. 
If a gamma ray line could be identified, it would constitute a "smoking 
gun" for dark matter annihilations. By definition, however, WIMPs do not 
annihilate through tree level processes to final states containing photons 
(if they did, they would be EMIMPs rather than WIMPs). On the other 
hand, they may be able to produce final states such as 77, 7Z or 7/1 through 
loop diagrams. Neutralinos, for example, can annihilate directly to 77 [57] 
or 7Z [58] through a variety of charged loops. These final states lead to 
gamma ray lines with energies of = TOdm and = mdm{^ — 'm%/'iiTi^jn)^ 
respectively. Such photons are produced in only a very small fraction of 
neutralino annihilations, however. The largest neutralino annihilation cross 
sections to 77 and are about 10~^® cm^/s, and even smaller values are 
more typical [59]. 

The Galactic Center has long been considered to be one of the most 
promising regions of the sky in which to search for gamma rays from dark 
matter annihilations [59, 60]. The prospects for this depend, however, on 
a number of factors including the nature of the WIMP, the distribution of 
dark matter in the region around the Galactic Center, and our ability to 
understand the astrophysical backgrounds present. 
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The gamma ray flux from dark matter annihilations is given by 

1 . . dN-y 1 

^7 47rr?^jf jios 



^^{E^,^) = -< axx\v\ > / p'(r)d/(^)#. (34) 



Here, < cxjc |w| > is the WIMP's annihilation cross section (times relative 
velocity), V is the angle observed relative to the direction of the Galac- 
tic Center, p(r) is the dark matter density as a function of distance to 
the Galactic Center, and the integral is performed over the line-of-sight. 
dN^/dE-y is the gamma ray spectrum generated per WIMP annihilation. 
Averaging over a solid angle centered around a direction, ip, we arrive at 

. .sxio-.-c.-^.-.g.( 3;;-i;i>^J (i^)^(A.,„.a 

(35) 

where AfJ is the solid angle observed. The quantity J(Af2, tp) depends only 
on the dark matter distribution, and is the average over the observed solid 
angle of the quantity 



J(V>) = -J— f I — -) [ p-^{r{l,'4,))dl. (36) 
8.5kpcV0.3GeV/cm3y Jio/ 



J{'ip) is normalized such that a completely flat halo proflle, with a density 
equal to the value at the solar circle, integrated along the line-of-sight to 
the Galactic Center would yield a value of one. In dark matter distribu- 
tions favored by N-body simulations, however, this value is much larger. A 
commonly used parameterization of halo profiles is given by 

where R ^ 20 kpc is the scale radius and po is fixed by imposing that the 
dark matter density at the distance of the Sun from the Galactic Center 
is equal to the value inferred by observations (-^0.3 GeV/cm^). Among 
the most frequently used parameterizations is the Navarro-Frenk- White 
(NFW) profile, which is described by a = 1, /3 = 3 and 7 = 1 [61]. When 
considering the region of the Galactic Center, the most important feature 
of the halo profile is the inner slope, 7. In some halo profiles, this slope can 
be considerably steeper than the value used in the NFW parameterization. 
For example, the Moore et al. profile is described by a = 1.5, (3 = 3, 
7 = 1.5 [62]. Note that for 7 > 1.5, the integral in Eq. (36) diverges. To 
avoid this behavior, one must truncate the profile within very small galactic 
radii. 
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The Narvarro-Frenk- White and Moore et al. profiles lead to values of 

J{AVt = 10~^ sr, ^ = Q) ^ 10^ and ^ 10^, respectively. And although 
these profiles serve as useful benchmarks, they certainly do not exhaust the 
range of possibilities. For a number of reasons, it is very difficult to predict 
the dark matter distribution in the highly important inner parsccs of the 
Galaxy. Firstly, the resolution of N-body simulations is limited to scales of 
approximately ~ 10^ parsecs or so. Furthermore, the gravitational potential 
in the inner region of the Milky Way is dominated not by dark matter, but 
by baryons, whose effects arc not generally included in such simulations. 
The precise impact of baryons on the dark matter distribution is difficult to 
predict, although an enhancement in the dark matter annihilation rate due 
to adiabatic compression is generally expected [63] . The adiabatic accretion 
of dark matter onto the central supermassive black hole may also lead to 
the formation of a density spike in the dark matter distribution. If present, 
such a spike would result in a very high dark matter annihilation rate [64] . 

The recent observation of a bright, very high energy gamma ray source 
in the region of the Galactic Center by HESS and other ground based At- 
mospheric Cerenkov Telescopes [65] has made efforts to identify gamma 
rays from dark matter annihilations more difficult. This source appears to 
be coincident with the dynamical center of the Milky Way (Sgr A*) and has 
no detectable angular extension (less than 1.2 arcminutes). Its spectrum is 
well described by a power-law, dNj/dE^ a E~", where a « 2.25 over the 
range of 160 GeV to 20 TeV, and thus appears to be inconsistent with a 
dark matter interpretation. Although this gamma ray source represents a 
formidable background for experiments searching for dark matter annihila- 
tion radiation [66] , it may be possible to reduce the impact of this and other 
backgrounds by studying the angular distribution of gamma rays from this 
region of the sky [67] . 

Telescopes potentially capable of detecting gamma rays from dark mat- 
ter annihilations include the satellite based Fermi gamma ray space tele- 
scope (formerly known as GLAST), and a number ground based Atmo- 
spheric Cerenkov Telescopes, including HESS, MAGIC and VERITAS. 
These two classes of experiments play complementary roles in the search 
for dark matter. On one hand, Fermi will continuously observe a large 
fraction of the sky, but with an effective area far smaller than possessed by 
ground based telescopes. Ground based telescopes, in contrast, study the 
emission from a small angular field, but with far greater exposure. Fur- 
thermore, while ground based telescopes can only study gamma rays with 
energy greater than ~100 GeV, Fermi will be able to directly study gamma 
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Figure 7. The projected exclusion limits at the 95% confidence level from five years 
of observation with the Fermi gamma ray space telescope (formerly known as GLAST) 
on the WIMP annihilation cross section as a function of its mass, for the case of an 
NFW halo profile. The region above the dotted line is already excluded by EGRET [68] . 
The dashed and solid lines show the projections for Fermi for an assumed isotropic 
diffuse background and in the limiting case in which the astrophysical background has 
the same angular distribution as the dark matter signal, respectively. Also shown are 
points representing a random scan of supersymmetric models. Figure from Ref. [67]. 

rays with energies over the range of 100 IVIeV to 300 GeV. 

In Fig. 7, the sensitivity of Fermi to dark matter annihilations in the 
Galactic Center region is shown for the case of an NFW halo profile and 
a WIMP annihilating to W^W^ . For this halo profile, WIMPs with an 
annihilation cross sections of < CTxjflwl >~ 3 x 10~^^ cm'^/s are near the 
threshold for detection by Fermi. For WIMPs with approximately this cross 
section, halo profiles more cuspy than NFW are thus likely to be observable 
by Fermi, whereas less dense profiles are unlikely to lead to an identifiable 
signal. 

If the dark matter density in the inner parsecs of the Milky Way is not 
particularly high, or if the astrophysical backgrounds turn out to be par- 
ticularly foreboding, the prospects for identifying dark matter annihilation 
radiation from the Galactic Center may be quite unfavorable. In this case, 
regions of the sky away from the Galactic Center may be more advanta- 
geous for dark matter searches. In particular, dwarf spheroidal galaxies 
within and near the Milky Way provide an opportunity to search for dark 
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matter annihilation radiation with considerably less contamination from 

astrophysical backgrounds. The flux of gamma rays from dark matter an- 
nihilations in such objects, however, is also expected to be lower than from 
a cusp in the center of the Milky Way [69]. As a result, planned experi- 
ments are likely to observe dark matter annihilation radiation from dwarf 
galaxies only in the most favorable range of particle physics models. Alter- 
natively, Fermi may also be sensitive to dark matter annihilations taking 
place throughout the halo of the Milky Way, or throughout the cosmological 
distribution of dark matter [70]. 

5.2. Charged Cosmic Rays From WIMP Annihilations 

In addition to gamma rays, WIMP annihilations throughout the galac- 
tic halo are expected to create charged cosmic rays, including electrons, 
positrons, protons and antiprotons. Unlike gamma rays, which travel along 
straight lines, charged particles move under the influence of the Galactic 
Magnetic Field, diffusing and steadily losing energy, resulting in a diffuse 
spectrum at Earth. By studying the spectrum of these particles, it may be 
possible to identify signatures of dark matter annihilations. In fact, multiple 
experiments have recently announced results which have been interpreted 
as possible products of WIMPs. 

The PAMELA experiment, which began its three-year satellite mis- 
sion in June of 2006, recently reported an anomalous rise in the cosmic 
ray positron fraction (the positron to positron-plus-electron ratio) above 
10 GeV [71], confirming earlier indications from HEAT [72] and AMS- 
01 [73]. Additionally, the ATIC balloon experiment has recently published 
data revealing a feature in the cosmic ray electron (plus positron) spec- 
trum between approximately 300 and 800 GeV, peaking at around 600 
GeV [74]. These measurements are each shown in Fig. 8, compared to the 
standard astrophysical predictions. These observations suggest the pres- 
ence of a relatively local (within ^1 kpc) source or sources of energetic 
cosmic ray electrons and positrons. Furthermore, in addition to the ob- 
servations of PAMELA and ATIC, the WMAP experiment has revealed 
an excess of microwave emission from the central region of the Milky Way 
which has been interpreted as synchrotron emission from a population of 
electrons/positrons with a hard spectral index [75, 76]. Taken together, 
these observations suggest that energetic electrons and positrons are sur- 
prisingly ubiquitous throughout our galaxy. 

Although the origin of these electrons and positrons is not currently 
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Figure 8. The positron fraction above 10 GeV as measured by the PAMELA experiment 
(top) [71] and the electron-plus-positron spectrum as measured by ATIC (bottom) [74]. 
In each frame, the range of standard astrophysical expectations is shown for compari- 
son. These observations clearly require an additional source of energetic electrons and 
positrons. Figures generated by Melanie Simet. 
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known, interpretations of the observations have focused two possibihties: 

emission from pulsars [77], and dark matter annihilations [78, 79]. In or- 
der for dark matter annihilations throughout in the Milky Way halo to 
produce a spectrum with a shape similar to that observed by PAMELA 
and ATIC, however, a large fraction of the annihilations must proceed to 
electron-positron pairs, or possibly to or t'^t^ [78]. Furthermore, 

WIMPs annihilating to other final states typically exceed the observed flux 
of cosmic ray antiprotons if normalized to generate the PAMELA and ATIC 
signals [80]. 

Dark matter particles which annihilate directly to e^e~ are predicted 
to generate a distinctive feature in the cosmic ray electron spectrum: an 
edge that drops off suddenly at = mx- In contrast, pulsars and other 
astrophysical sources of cosmic ray electrons arc expected to produce spec- 
tra which fall off more gradually. Although the current data from ATIC is 
not detailed enough to discriminate between a feature with a sudden edge 
(dark matter-like) or graduate cutoff (pulsar-like), such a discrimination 
could become possible if the electron spectrum were measured with greater 
precision. Interestingly, such a measurement should be possible for ground 
based gamma ray telescopes such as HESS or VERITAS [81]. 

Once electrons and positrons arc injected into the local halo through 
dark matter annihilations (or from pulsars), they propagate under the in- 
fluence of the Galactic Magnetic Field, gradually losing energy through 
synchrotron emission and through inverse Compton scattering with radia- 
tion fields. At energies of a few GcV and higher, the resulting spectrum at 
Earth can be calculated by solving the diffusion- loss equation [82]: 



d dUe 
dtdK 



+ 



U TP -\ 



+ Q{E„x), (38) 



where dn^/dE^ is the number density of positrons/electrons per unit energy, 
K{Ee,x) is the diffusion constant, b{Ee,x) is the rate of energy loss, and 
Q{Ef,, x) is the source term, which contains all of the information regarding 
the dark matter annihilation modes, cross section, and spatial distribution. 
As we expect the cosmic ray distribution from dark matter annihilations 
to be in or near the steady state limit, the left hand side of the equation is 
generally set to zero. The energy dependance of the diffusion constant is 
typically parameterized by K(Ep) = KqE^- 

In the relativistic limit, the energy loss rate resulting from inverse Comp- 
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ton scattering and synchrotron emission is given by 

« 1.02 X lO-i^^GeV/s ( ^£2±±^^ x ( , (39) 
\ eV/cm'^ J \GeVy 

where ctt is the Thompson cross section and prad and pmag arc the energy 
densities of radiation (including starlight, cosmic microwave background 
and emission from dust) and the galactic magnetic field, respectively. 

To find a solution to the diffusion-loss equation, a set of boundary con- 
ditions must also be adopted. For example, under the assumption of cylin- 
drical symmetry, the distance from the galactic plane at which electrons 
and positrons freely escape, L, must be selected. 

The diffusion parameters Kq, a and L can be constrained by studying 
the spectra of various species of cosmic ray nuclei. In particular, by studying 
the ratios of cosmic ray secondaries-to-primaries (most importantly boron- 
to-carbon) and unstable primarics-to-stablc secondaries (Be^'^-to-Be^), in- 
formation can be inferred regarding the size of the diffusion zone and the 
length of time that cosmic rays are confined to the galaxy. Current cosmic 
ray data favor values of approximately Kq ~ 10^* cm^/s, a ~ 0.4 — 0.5 and 
i ~ 1 — 10 kpc. For an excellent review of this subject, sec Ref. [83]. 

The spectral shape of the source term, Q, depends on the leading annihi- 
lation modes of the WIMP in the low velocity limit. Neutralinos and other 
WIMPs which annihilate primarily to combinations of heavy fcrmions and 
Higgs or gauge bosons generally produce somewhat soft spectra. In con- 
trast, Kaluza-Klein dark matter and other WIMPs which annihilate sig- 
nificantly to electrons and muons are predicted to generate a considerably 
harder spectrum [84]. 

The normalization of the source term is determined by the annihilation 
cross section of the WIMP and the spatial distribution of the dark matter 
in the galactic halo. In particular, small scale inhomogeneities in the dark 
matter distribution can boost the average annihilation rate. The factor by 
which inhomogeneities enhance the resulting cosmic ray flux is called the 
"boost factor" , and is defined as 

boost factor = / ^ (40) 

where the integral is performed over the volume that contributes to the 
observed spectrum. The results of N-body simulations lead us to expect 
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that this quantity could be as large as 5 to 10. Although considerably larger 
boost factors arc not impossible, they would be somewhat surprising. 

In order to generate the signals observed by PAMELA and ATIC with 
a WIMP annihilating with an annihilation cross section of crw « 3 x 10~^^ 
cm^/s (as required to thermally produce the observed dark matter abun- 
dance with a fixed value of rrv) an annihilation rate many hundred times 
larger than is naively expected would be required. This could potentially be 
accommodated in several ways. Firstly, the boost factor may be larger than 
expected. Secondly, the dark matter particles may have been generated 
through a non-thermal mechanism, allowing for considerably larger annihi- 
lation cross sections. Thirdly, WIMPs interacting through the exchange of 
very light particles can annihilate through non-perturbative processes such 
that av oc v~^. In this case, the annihilation cross section in the present 
galaxy (with v ~ 10""^ c) would be much larger than the corresponding 
cross section at the time of thermal freeze-out (when v ~ 10~^ c) [85]. 

5.3. Neutrinos From WIMP Annihilations in the Sun 

As the Solar System moves through the halo of the Milky Way, WIMPs 
become swept up by the Sun. Although dark matter particles interact only 
weakly, they occasionally scatter elastically with nuclei in the Sun and lose 
enough momentum to become gravitationally bound. Over the lifetime of 
the Sun, a sufHcient density of WIMPs can accumulate in its center so 
that an equilibrium is established between their capture and annihilation 
rates. The annihilation products of these WIMPs include neutrinos, which 
escape the Sun with minimal absorption, and thus potentially constitute an 
indirect signature of dark matter. Such neutrinos can be generated through 
the decays of heavy quarks, gauge bosons, and other products of WIMP 
annihilation, and then proceed to travel to Earth where can be efficiently 
identified using large volume neutrino detectors. 

Beginning with a simple estimate, we expect WIMPs to be captured in 
the Sun at a rate approximately given by: 

^<l>x{MQ/mp)axp, (41) 

where (j)x is the flux of WIMPs in the Solar System, Mq is the mass of the 

Sun, and axp is the WIMP-proton elastic scattering cross section. Rea- 
sonable estimates of the local distribution of WIMPs leads to a capture 
rate of C® ~ lO^Osec"! x (lOOGeV/mx) (o-xp/lQ-^pb). This neglects, 
however, a number of potentially important factors, including the gravita- 
tional focusing of the WIMP flux toward the Sun, and the fact that not 
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every scattered WIMP will ultimately be captured. Taking these and other 
effects into account leads us to a solar capture rate of [86]: 



C- . 1.3 X 10- sec- ( 3 ) ( '-^^^ 

V0.3GcV/cm'^/ V wiocai , 

100GcV \ ^ / A (axi,SD+ (Txi,si) 
mx J ^\ 10-6 pb j 



(42) 



where piocai is the local dark matter density and ■uiocai is the local rms veloc- 
ity of halo dark matter particles, ctxi.sd and ctxi.si are the spin-dependent 
and spin-independent elastic scattering cross sections of the WIMP with 
nuclei species i, and Ai is a factor denoting the relative abundance and 
form factor for each species. In the case of the Sun, Ah « 1.0, Aae ~ 0.07, 
and Ao ~ 0.0005. The quantity S contains dynamical information and is 
given by 



S{x) 

where 



f(x)3/2 

1 + F(a:)3/2 



2/3 

(43) 



2 (X - 1)^ V'^local/ 

and Vesc ~ 1156 km/s is the escape velocity of the Sun. Notice that for 
WIMPs much heavier than the target nuclei 5 oc 1 /mx, leading the capture 
rate to be suppressed by two factors of the WIMP mass. In this case 
{mx ^ 30 GeV), we can write the capture rate as: 

C- . 3.35 X 10- sec- ( ^^-] \ ('J^] ' ('-^] ' 
V0.3GeV/cm3y \ wiocai J \ mx 

f ajfH,SD+ cTjfH,si + 0.07 axHe,si + 0.0005 S'(^)axo,si 
^ 10-6 pb 

If the capture rate and annihilation cross sections are sufficiently large, 
equilibrium will be reached between these processes. For N{t) WIMPs in 
the Sun, the rate of change of this quantity is given by 

N{t) = C® - A®N{tf - E^N, (46) 

where C® is the capture rate described above and A® is the annihilation 
cross section times the relative WIMP velocity per volume. E® is the 
inverse time for a WIMP to escape the Sun via evaporation. Evaporation 
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is highly suppressed for WIMPs heavier than a few GeV [87, 88]. ^4® can 
be approximated by 

A- = M, (47) 

where Kff is the effective volume of the core of the Sun determined roughly 
by matching the core temperature with the gravitational potential energy 
of a single WIMP at the core radius and is given by [87, 88] 

K. = 5.7xlO-cm3fi^V'\ (48) 

Neglecting evaporation, the present WIMP annihilation rate is given by 

r = ^A'^Nitef = ^ C® tanh^ (\/C®A® t©) , (49) 

where tQ « 4.5 billion years is the age of the solar system. The annihilation 
rate is maximized when it reaches equilibrium with the capture rate. This 
occurs when 

\/C©A®iQ > 1 . (50) 

If this condition is met, the final annihilation rate is determined entirely 
by the capture rate and has no further dependence on the dark matter 
particle's annihilation cross section. 

Through their annihilations, WIMPs can generate neutrinos through 
a variety of channels. Annihilations to heavy quarks, tau Icptons, gauge 
bosons and/or Higgs bosons can each generate energetic neutrinos in their 
subsequent decays [89]. In some models, WIMPs can also annihilate directly 
to neutrino- antineutrino pairs. Annihilations to light quarks or muons, 
however, do not contribute to the high energy neutrino spectrum, as these 
particles come to rest in the solar medium before decaying. 

Neglecting the effects of oscillations and interactions with the solar 
medium, the spectrum of neutrinos from WIMP annihilations to a final 
state, YY, is given by: 

dN, 1 /•^-/'^(i-W 1 dE' fdN, , 

(51) 



dE^ 2 J E,j^(i+I3) E' \dE„ J YY 



where 7 = nix/mY, P = \/l — 7"^, and {dN^/dE^)YY is the spectrum of 
neutrinos produced in the decay of an Y at rest. 

Gauge bosons produced in WIMP annihilations produce the most ener- 
getic neutrinos through their decays, W ^ If, Z ^ vv, but also produce 
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neutrinos through the subsequent decays of muons, taus and other parti- 
cles. Tau Icptons produce neutrinos through a variety of channels, including 
through the semi-leptonic decays r fivv, evv, and the hadronic decays 
r TTU, Ku, TTTTV, and ■k-kttv. Top quarks decay to a and a bottom 
quark essentially 100% of the time, each of which can generate neutrinos 
in their subsequent decay. For bottom and charm quarks, only the semi- 
leptonic decays contribute to the neutrino spectrum (with the exception 
of the neutrinos resulting from taus and c-quarks produced in decays of 
6-quarks). In the case of h and c quark decays, the process of hadronization 
reduces the fraction of energy that is transferred to the resulting neutrinos 
and other decay products. 

Once produced in the Sun's core, neutrinos then propagate through the 
solar medium and to the Earth. Over this journey, they can potentially be 
absorbed, lose energy, and/or change flavor. In particular, charged current 
interactions of electron and muon neutrinos in the Sun lead to their ab- 
sorption. The probability of absorption taking place can be estimated by 
1 — exp(— £^y/£abs)j where i^abs is approximately 130 GeV for electron or 
muon neutrinos and 200 GeV for electron or muon antineutrinos. Absorp- 
tion, therefore, only plays an important role for relatively heavy WIMPs. 

The effect of charged current interactions on tau neutrinos in the Sun is 
somewhat more complicated. The tau leptons produced in such interactions 
quickly decay and thus regenerate the absorbed tau neutrino, albeit with 
a reduced energy. Neutral current interactions of all three neutrino flavors 
similarly reduce the neutrinos' energy without depleting their number. 

Vacuum oscillations lead to the full mixing of muon and tau neutrinos 
over their propagation to the Earth, making the observed muon neutrino 
spectrum (which is the flux relevant for detection in neutrino telescopes) 
effectively the average of the muon and tau flavors prior to mixing. Electron 
neutrinos can also oscillate into muon flavor through matter eff'ects in the 
Sun (the MSW effect). Electron antineutrinos can generally be neglected, 
as their oscillations to muon or tau flavors are highly suppressed [90]. 

Program such as DarkSUSY [29], which include effects such as 
hadronization, absorption, regeneration, and oscillations, are very useful 
in making detailed predictions for the neutrino spectrum resulting from 
WIMP annihilations in the Sun. 

Once they reach Earth, neutrinos can potentially be detected in large 
volume neutrino telescopes. Neglecting oscillations and solar absorption, 
the muon neutrino spectrum at the Earth resulting from WIMP annihila- 
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tions in the Sun is given by 



(52) 



where Cq is the WIMP capture rate in the Sun, F^q is the non-equihbrium 
suppression factor (i^Eq — 1 for capture-annihilation cquihbrium), Des is 
the Earth-Sun distance and (^^^)^"'' is the muon neutrino spectrum from 
the Sun per WIMP annihilating. 

Muon neutrinos can produce muons through charged current interac- 
tions with ice or water nuclei inside or near the detector volume of a high 
energy neutrino telescope. The rate of neutrino-induced muons observed in 
a high energy neutrino telescope is estimated by 

iVevents ^ J J ^ ^ [E,, , v) [R^^{E^,) + L] A,s dE,^ dy 

+ J i^^. ' y) + L] Aff dE^^ dy, (53) 

where Oi, (cxp) is the neutrino- nucleon (antineutrino-nucleon) charged cur- 
rent interaction cross section, (1 — y) is the fraction of neutrino/antineutrino 
energy which goes into the muon, is the effective area of the detector, 
R^{E^) is the distance a muon of energy, (1 — y) E,^, travels before falling 
below the energy threshold of the experiment (ranging from approximately 
1 to 100 GeV), called the muon range, and L is the depth of the detector 
volume. The muon range in water/ice is approximately given by 

2.0 0.0042 S^(GcV) 



R^{E^) K. 2.4km X In 



2.0 0.0042 i;t'i'-(GeV) 



(54) 



where , E^^"^ is the threshold of the experiment. 



When completed, the IceCube experiment will possess a full square 
kilometer of effective area and kilometer depth, and will be sensitive to 
muons above approximately 50 GcV [91]. The Deep Core extension of 
Icecube will be sensitive down to 10 GeV. The Super-Kamiokande detector, 
in contrast, has 10~^ times the effective area of IceCube and a depth of 
only 36.2 meters [92]. For low mass WIMPs, however, Super-Kamiokande 
benefits over large volume detector such as IceCube by being sensitive to 
muons with as little energy as ~1 GeV. 

The spectrum and flux of neutrinos generated in WIMP annihilations 
depends on the annihilation modes which dominate, and thus arc model 
dependent. As long as the majority of annihilations proceed to final states 
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Figure 9. The event rate in a kilometer-scale neutrino tcleseope such as IceCube as a 
function of the WIMP's cfTcctive elastic scattering cross section in the Sun for a variety 
of annihilation modes. The effective elastic scattering cross section is defined as (Joff = 
o-XH,SD + crxH,Sl + 0.07crxHe,Sl + 0.0005 5(mx /mo) crxo, SI (see Eq. 45). The dashed, 
solid and dotted lines correspond to WlMPs of mass 100, 300 and 1000 GeV, respectively. 
A 50 GeV muon energy threshold and an annihilation cross section of 3 x 10~^^ cm^ 
have been adopted. 



such as hb, tt, r+r~, ZZ, or some combination of Higgs and gauge 

bosons, however, the variation between different final states is not dramatic. 
In Fig. 9, we plot the event rate in a kilometer-scale neutrino telescope such 
as IceCube as a function of the WIMP's effective elastic scattering cross 
section for four possible annihilation modes. The effective elastic scattering 
cross section used here is defined as a^g = ctxh.sd + fxH.si + 0-07crxHe,Si + 
0.0005 5(TOx/mo)axo,si (see Eq. 45). 

The elastic scattering cross section of a WIMP is constrained by the 
absence of a positive signal in direct detection experiments. As described in 
Sec. 4, the strongest limits on the WIMP-nucleon spin-independent elastic 



39 



scattering cross section have been made by the CDMS [42] and XENON [43] 
experiments. These results cxchide spin-independent cross sections larger 
than approximately 5 x 10"^ pb for a 25-100 GeV WIMP or 2 x 10"^ pb 
X {mx /500 GeV) for a heavier WIMP. 

With these results in mind, consider as an example a 300 GeV WIMP 
with an elastic scattering cross section with nuclcons which is largely spin- 
independent. With a cross section near the CDMS bound, say 1 x 10~^ 
pb, we can determine from Fig. 9 the corresponding rates in a kilometer- 
scale neutrino telescope, such as IceCube. Sadly, we find that this cross 
section yields less than 1 event per year for annihilations to 66, about 2 
events per year for annihilations to W~^W~ or tt and about 8 per year 
for annihilations to r+r", none of which are likely to be sufficient to be 
identified over the background of atmospheric neutrinos by IceCube or other 
planned experiments. Clearly, WIMPs that scatter with nucleons mostly 
through spin-independent interactions are not likely to be detected with 
IceCube or other planned neutrino telesc;opc;s. 

The same conclusion is not reached for the case of spin-dependent 
scattering, however. The strongest bounds on the WIMP-proton spin- 
dependent cross section have been placed by the COUPP [49] and KIMS [51] 
collaborations. These constraints are approximately 7 orders of magnitude 
less stringent than those corresponding to spin- independent couplings, how- 
ever. As a result, a WIMP with a largely spin-dependent scattering cross 
section with protons may be capable of generating large event rates in high 
energy neutrino telescopes. Again considering a 300 GeV WIMP with a 
cross section near the experimental limit. Fig. 9 suggests that rates as high 
as ~10^ per year could be generated if purely spin-dependent scattering 
contributes to the capture rate of WIMPs in the Sun. 

Spin-dependent, axial-vector, scattering of neutralinos with the quarks 
or gluons within a nucleon is made possible through the t-channel exchange 
of a Z-boson, or the s-channel exchange of a squark. Although the cross 
sections for these processes can vary dramatically depending on the details 
of the supersymmetric model under consideration, it is often the case that 
the neutralino's spin-dependent cross section is considerably larger than 
its corresponding spin-independent interaction. In particular, very large 
spin-dependent cross sections (ctsd ^ 10~^pb) are possible even in models 
with very small spin-independent scattering rates. Such a model would go 
easily undetected in all planned direct detection experiments, while still 
generating on the order of ^ 1000 events per year at IceCube. 

In Fig. 10, we demonstrate this by plotting the rate in a kilometer-scale 
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Figure 10. The rate of events at a kilometer-scale neutrino telescope such as IceCube 
from neutralino annihilations in the Sun, as a function of the neutralino's spin-dependent 
elastic scattering cross section with protons. Each point shown is beyond the reach of 
current direct detection experiments. 

neutrino telescope from WIMP annihilation in the Sun verses the WIMP's 
spin-dependent cross section with protons. In this figure, each point shown 
is beyond the reach of present direct detection experiments. We thus con- 
clude that neutralinos may be observable by IceCube while remaining un- 
observed by current and near future direct detection experiments. 

A neutralino which has a large spin-dependent cross section generally 
has a sizable coupling to the Z, and thus has a large higgsino component. 
In particular, the spin-dependent scattering cross section through the ex- 
change of a Z is proportional to the square of the quantity |iVi3p — |A''i4p. 
As a result, neutralinos with a few percent higgsino fraction or more are 
likely to be within the reach of IceCube [93, 94]. This makes the focus 
point region of supersymmetric parameter space especially promising. In 
this region, the lightest neutralino is typically a strong mixture of bino and 
higgsino components, often leading to the prediction of hundreds of events 
per year at IceCube. 

Considering Kaluza-Klein dark matter, the range of elastic scattering 
cross sections predicted are quite challenging to reach with direct detec- 
tion experiments, but are more favorable for detection using neutrino tele- 
scopes. The spin-independent S^^-'-nucleon cross section, which is gener- 
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ated through the exchange of KK quarks and the Higgs boson, is rather 

small and well beyond the sensitivity of current or upcoming direct detec- 
tion experiments. The spin-dependent scattering cross section for the B^^^ 
with a proton, however, is considerably larger and is given by [54] 

800GeV\^ /O.l 



where Ag = (m^ci) — mg(i))/mgii) is fractional shift of the KK quark 
masses over the -B^^^ mass, which is expected to be on the order of 10%. 
The /S"s parameterize the fraction of spin carried by each variety of quark 
within the proton. 

In addition to this somewhat large spin-dependent scattering cross sec- 
tion, the annihilation products of the B^^^ are very favorable for the pur- 
poses of generating observable neutrinos. In sharp contrast to neutralinos, 
approximately 60% of B^^^ annihilations generate a pair of charged leptons 
(20% to each type). Although most the remaining 40% of annihilations 
produce up-type quarks, about 4% generate neutrino pairs directly. The 
neutrino and tau lepton final states each contribute substantially to the 
event rate in a neutrino telescope. 

Taken together, this leads to the prediction of a fairly high neutrino 
event rate from Kaluza-Klein dark matter annihilating in the Sun. In par- 
ticular, for masses in the 500-1000 GeV range and a 10% mass degeneracy 
with the Kaluza-Klein quarks, we expect ^10-1000 events per year in a 
kilometer-scale neutrino telescope such as IceCube [95]. 

Currently, the strongest constraints on the neutrino flux from WIMPs 
annihilating in the Sun have been placed by the IceCube (for nix ^ 200 
GeV) and Super-Kamiokande {mx ^ 200 GeV) collaborations. The current 
IceCube limit constrains the neutrino- induced muon rate from dark matter 
to be less than ~400-500 per square kilometer per year. Ultimately, IceCube 
is expected to reach a sensitivity about an order of magnitude below this 
current level. 



6. Signals, Hints, and... Otherwise 

Over the past several years, a number of observations have been interpreted 
as possible products of dark matter annihilations. In this section, I take 
the opportunity to summarize and discuss some of these observations. 
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6.1. The PAMELA and ATIC Excesses 

As I have already discussed the recent observations of the ATIC [74] and 
PAMELA [71] experiments in Sec. 5.2, in this section I will only briefly 
summarize some of the arguments in favor of and against these signals 
being likely products of dark matter annihilations. 

Although the observations of ATIC and PAMELA support a very com- 
pelling case that a powerful source of energetic positrons and electrons is 
present within ~1 kpc of the Solar System, the nature of this source or 
sources is not yet clear. In order to produce the PAMELA and ATIC sig- 
nals through the annihilations of dark matter distributed throughout the 
galactic halo, the WIMPs must annihilate dominantly to charged leptons 
(to produce a sufficiently hard spectrum, and to avoid the overproduction of 
cosmic ray antiprotons) [78, 79, 80]. Furthermore, a very large annihilation 
rate is also required - hundreds or thousands of times higher than is naively 
expected for a thermal relic. To accommodate such a high annihilation rate, 
we must require that either the WIMPs possess a very large annihilation 
cross section (which, in turn, requires a non-thermal mechanism for their 
production in the early universe), or that local inhomogeneities in the dark 
matter distribution boost the annihilation rate far more efficiently than N- 
body simulations would lead us to expect. Large annihilation cross sections 
might also result from non-pcrturbativc processes in some models [85]. Al- 
ternatively, in the relatively unlikely event that a large and dense clump of 
dark matter happened to reside within ~1 kpc of the Solar System, a suf- 
ficiently high annihilation rate, hard electron/positron spectrum, and low 
antiproton flux could plausibly be generated [96]. 

Taken together, these considerations lead me to conclude that although 
the ATIC and PAMELA signals could potentially be explained by dark 
matter annihilations, such a scenario would require WIMPs which possess 
rather special properties, or that are distributed in a somewhat unlikely 
way. The leading astrophysical alternative for the origin of these signals is 
a nearby pulsar (or pulsars) [77]. Although pulsars are known to be sites 
of electron-positron pair production, it is not possible to reliably predict 
the spectrum of or total power injected from these objects. To accom- 
modate the PAMELA and ATIC observations, a nearby (within -^1 kpc) 
and somewhat young (10^-10^ years) pulsar or pulsars must have deposited 
a few percent or more of their total energy output in the form of a very 
hard spectrum of electron-positron pairs. While this appears to be a larger 
fraction of the total energy than was generally expected, it is certainly not 
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implausible. 

In the relatively near future, measurements by the Fermi gamma ray 
space telescopes, as well as ground based telescopes [81], should clarify this 
situation considerably. The origin of the PAMELA and ATIC signals will 
likely not remain a mystery for long. 

6.2. The WMAP Haze 

In addition to its measurements of the cosmic microwave background, 
data from the Wilkinson Microwave Anisotropy Probe (WMAP) has been 
used to provide the best measurements to date of the standard interstellar 
medium emission mechanisms, including thermal dust, spinning dust, ion- 
ized gas, and synchrotron. In addition to these expected foregrounds, the 
observations have revealed an excess of microwave emission in the inner 20° 
around the center of the Milky Way, distributed with approximate radial 
symmetry. This excess is known as the "WMAP Haze" [75]. 

Although the WMAP Haze was initially thought likely to be thermal 
bremsstrahlung (free-free emission) from hot gas (10'*K ^ T ^ 10^ K), 
this interpretation can be ruled out by the relative absence of the Ha re- 
combination line and X-ray emission. Other possible origins for this signal, 
such as thermal dust, spinning dust, and Galactic synchrotron as traced 
by low-frequency surveys, also seem unlikely. Alternatively, it has been 
suggested that the WMAP Haze could be a product of dark matter an- 
nihilations [76]. In particular, annihilating dark matter particles produce 
relativistic electrons and positrons which travel under the influence of the 
Galactic magnetic field. As they do, they will emit synchrotron photons, 
which naturally fall within the frequency range measured by WMAP. 

The angular distribution of the Haze can be used to constrain the shape 
of the required dark matter halo profile. In particular, the morphology of 
the Haze is consistent with originating from dark matter distributed as 
p{r) oc r^'^-'^ within the inner kiloparsecs of our galaxy [76]. This slope 
falls between those predicted by the NFW, p(r) oc r~^, and Moore et al., 
p{r) (X r~^-^, halo profiles. The annihilation cross section required of a 
~100-1000 GeV WIMP to produce the observed intensity of the WMAP 
Haze is of the same order magnitude as the value predicted for a simple 
thermal relic. No large boost factors are needed to generate this signal. 

It is also interesting to note that the dark matter halo profile and annihi- 
lation cross section required to generate the WMAP Haze with dark matter 
imply a flux of prompt gamma rays from the Galactic Center region that 
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is within the reach of the upcoming Fermi gamma ray space telescope [97]. 

Additionally, the upcoming Planck satellite will provide substantially im- 
proved measurements of the spectrum and angular distribution of the Haze. 



6.3. DAMA 's Annual Modulation 

The direct detection experiment DAMA has reported evidence for an an- 
nual modulation in its rate of nuclear recoil events [48] . It has been claimed 
that this signal is the result WIMP interactions, which result from varia- 
tions in the relative velocity of the Earth with respect to the dark matter 
halo as the Earth orbits the Sun. This effect is predicted to lead to a vari- 
ation in the flux of dark matter particles and their velocity distribution, 
with expected extrema occurring at June 2 and December 2. The DAMA 
experiment observes a maximum rate at low nuclear recoil energies on May 
24, plus or minus 8 days, and have accumulated enough data to put the 
significance of the observed modulation at approximately 8a. The collab- 
oration has not been able to identify any other systematic effects capable 
of producing this signal. The claim that this signal is the resiilt of dark 
matter interactions has been controversial, in part because a number of 
other experiments appear to be in direct conflict with the DAMA result. 

Several studies have attempted to reconcile the DAMA modulation sig- 
nal with the null results of other direct-detection experiments [98]. In 
particular, an elastically scattering WIMP with a mass in the several GeV 
range can satisfy the results of DAMA while remaining marginally consis- 
tent with the null resuhs of CDMS [42], CRESST [46], CoGeNT [47], and 
XENON [43]. The allowed parameter region depends crucially on the oc- 
currence of an effect known as channeling in the Nal crystals of the DAMA 
apparatus [99]. 

Another possibility to have been proposed is that the DAMA signal 
might arise from a WIMP which does not scatter with nuclei elastically, 
but instead scatters inelastically, leaving the interaction in the form of a 
slightly heavier state (with a mass splitting on the order of 100 keV) [100]. 
For kinematic reasons, this scenario allows for the efficient scattering of 
the WIMPs with iodine nuclei in DAMA, while suppressing the scattering 
rate off of germanium and other comparatively light nuclei used in other 
experiments. 
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6.4. The INTEGRAL 511 keV Line 

In 2003, the SPI spectrometer onboard the INTEGRAL sateUite confirmed 

the very bright emission of 511 kcV photons from the region of the Galac- 
tic Bulge, corresponding to an injection rate of approximately 3 x 10^^ 
positrons per second in the inner Galaxy [101]. This is orders of magnitude 
larger than the expected rate from pair creation via cosmic ray interactions 
with the interstellar medium. The signal appears to be approximately (but 
not perfectly [102]) spherically symmetric (with a full- width-half-maximum 
of approximately 6°), with little of the emission tracing the Galactic Disk. 
A stellar origin of this signal, such as type la supernovac, hypcrnovae or 
gamma ray bursts, would therefore also require a mechanism (such as sub- 
stantial coherent magnetic fields) by which the positrons could be trans- 
ported from the disk to throughout the volume of the Bulge [103]. Further- 
more, type la siipernovae do not inject enough positrons to generate the 
observed intensity of this signal [104]. It is possible, however, that hyper- 
novae [105] or gamma ray bursts [105, 106] might be capable of injecting 
positrons at a sufficient rate. A large population of several thousand X-ray 
binaries has also been proposed as a possible source of these photons [107]. 

Given the challenges involved with generating the observed 511 keV 
emission with astrophysical sources, it was suggested that this signal could 
potentially be the product of dark matter annihilations [108]. In order 
for dark matter particles to generate the observed spectral line width of 
this signal, however, their annihilations must inject positrons with energies 
below a few MeV [109]. This, in turn, implies that the dark matter's 
mass be near the 1-3 MeV range - much lighter than annihilating dark 
matter particles in most theoretically attractive models (for an interesting 
exception, see Ref. [110]). 

Although weakly interacting particles with masses smaller than a few 
GeV (but larger than ~1 MeV) tend to be overproduced in the early uni- 
verse relative to the measured dark matter abundance [16], this can be 
avoided if a new light mediator is introduced which makes dark matter an- 
nihilations more efficient [111]. For example, although neutralinos within 
the MSSM are required by relic abundance considerations to be heavier than 
-^20 GeV [112], they can be much lighter in extended supersymmetric mod- 
els in which light Higgs bosons can mediate neutralino annihilations [113]. 

For dark matter particles with MeV-scale masses to generate the mea- 
sured dark matter abundance, they must annihilate during the freeze-out 
epoch with a cross section of ctw ^ 3 x 10~^^ cm^/s. To inject the flux of 
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positrons needed to generate the signal observed by SPI/INTEGRAL, how- 
ever, an annihilation cross section fonr to five orders of magnitude smaller 
is required. Together, these requirements force us to consider dark matter 
particles with annihilation cross sections of the form, <jv ocv^. Such behav- 
ior can be found, for example, in the case of fermionic or scalar dark matter 
particles annihilating through a vector mediator. For such a dark matter 
particle to not be overabundant today, the mediating boson also must be 
quite light [114]. 

As an alternative explanation for the 511 keV line, it has been proposed 
that ~ 500 GeV dark matter particles could be coUisionally excited to states 
which are -^1 MeV heavier, which produce electron-positron pairs in their 
subsequent de-excitations to the ground state [115]. 

6.5. egret's Diffuse Galactic Spectrum 

The satellite-based gamma ray detector, EGRET, has measured the diffuse 

spectnim of gamma rays over the entire sky. When compared to conven- 
tional galactic models, these measurements appear to contain an excess at 
energies above approximately 1 GeV. This has been interpreted as evidence 
for dark matter annihilations in the halo of the Milky Way [116]. 

Among the most intriguing features of the observed EGRET excess is 
its similar spectral shape over all regions of the sky. Furthermore, this 
shape is consistent with that predicted from the annihilations of a 50-100 
GeV WIMP. There arc, however, some substantial challenges involved with 
interpreting the EGRET excess as a product of dark matter annihilations. 
In particular, to accommodate the required normalization for the annihila- 
tion rate in various regions of the Galaxy, the distribution of dark matter 
has to depart substantially from the predictions of standard dark matter 
halo profiles. In particular, Ref. [116] adopts a distribution which includes 
two very massive 10^° Mq) toroidal rings of dark matter near or within 
the Galactic Plane, at distances of approximately 4 and 14 kiloparsecs from 
the Galactic Center. The authors motivate the presence of these rings by 
observed features in the Galactic rotation curve, and suggest that they 
may be remnants of very massive dwarf galaxies which have been tidally 
disrupted. 

The other difficulty involved with the interpretation of the EGRET 
excess as dark matter annihilation radiation is the large flixx of antiprotons 
which is expected to be generated in such a scenario [117] . In particular, the 
flux of cosmic antiprotons produced is expected to exceed the measured flux 
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by more than an order of magnitude. To avoid this conclusion, one is forced 
to consider significant departures from standard galactic diffusion models. 
In particular, an anisotropic diffusion model featuring strong convection 
away from the Galactic Disk and a large degree of inhomogeneities in the 
local environment could reduce the cosmic antiproton flux to acceptable 
levels [118]. 

The dark matter interpretation of EGRET's measurement of the galac- 
tic diffuse spectrum has also been challenged on the grounds that the data 
could plausibly bo explained without the addition of an exotic component, 
from dark matter or otherwise. In particular, uncertainties in the cosmic 
ray propagation and diffusion model lead to considerable variations in the 
predicted diffuse gamma ray backgrounds [119]. More recently, it has also 
been suggested that the observed excess could also be the result of system- 
atic errors in EGRET's calibration [120]. 

The Fermi gamma ray space telescope should be able to quickly clarify 
the origin of the EGRET diffuse emission. 

6.6. EGRET's Diffuse Extragalactic Spectrum 

In addition to galactic emission, it has also been proposed that dark matter 
annihilation radiation might constitute a significant fraction of the extra- 
galactic (isotropic) diffuse gamma ray flux as measured by EGRET [121]. 

The origin of EGRET's diffuse extragalactic gamma ray background is 
currently unknown. Although all or most of the observed spectrum could 
be the product of astrophysical source such as blazars, much of the flux 
observed in the 1-20 GeV range could also plausibly be the result of dark 
matter annihilations taking place throughout the universe [70] . In particu- 
lar, the observed spectrum fits reasonably well the predictions for a WIMP 
with a mass of roughly 500 GeV. 

The dark matter annihilation rate needed to normalize to the diffuse 
flux measured by EGRET is, however, quite high and requires cither a very 
large dark matter annihilation cross section or dark matter halos which 
are very cuspy. In particular, if an NFW profile [61] is adopted for all 
halos throughout the universe, then an annihilation cross section 10^ to 
10'' times larger than is predicted for a thermal relic is required to generate 
the observed gamma ray fiux. If the dark matter distribution in the Milky 
Way is similar to that found in halos throughout the universe, however, 
then the gamma ray flux from the center of our galaxy would far exceed 
that which is observed [122]. To generate the isotropic diffuse flux observed 
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by EGRET without conflicting with observations of the Galactic Center, 

therefore, requires extremely cusped halo profiles in most or at least many 
of the galaxies throughout the universe, and a far less dense cusp in our 
own Milky Way. 

7. Summary and Outlook 

The next few years will be a very exciting time for research in particle 
dark matter. As direct detection experiments move closer to the ton-scale, 

many of the most attractive models of dark matter will become within 
their reach. The lack of a detection in these experiments would be capable 
of severely constraining the nature of supersymmetry or other TeV-scale 
physics containing a WIMP candidate. A confirmed detection of a WIMP 
would open the door to the age of precision dark matter studies, in which 
measurements of the WIMP's mass, interactions, and distribution would 
begin to be made. 

Indirect detection efforts are currently developing very rapidly. In par- 
ticular, much of what I have written here about electron, positron and 
gamma ray signals for dark matter will be hopelessly out of date even a 
year from now. But rapid progress makes for exciting times! With new 
data from the Fermi gamma ray space telescope on the horizon, the near 
future will likely hold many exciting results for indirect detection. 

The single most remarkable aspect of these lecture notes is that I have 
managed to write almost 50 pages about TeV-scale dark matter without 
mentioning the Large Hadron Collider. This should not be taken as an 
indication that the LHC is not important in the hunt for dark matter's 
identity. In contrast, I fully expect the next twenty years of particle physics 
(including particle-astrophysics) to be largely defined by what this incred- 
ible machine reveals to us. And while I remain largely agnostic regarding 
what the LHC is likely to discover, it will almost certainly give us insights 
into the nature of dark matter (even if that insight is that dark matter is 
not made up of WIMPs). 

Colliders and astrophysical experiments each provide very different and 
complementary types of information regarding dark matter. It is very un- 
likely that any single experiment or class of experiments will be sufficient 
to conclusively identify the particle nature of dark matter. The direct or 
indirect detection of the dark matter particles making up our galaxy's halo 
will probably not be able to provide enough information to reveal the un- 
derlying physics (supersymmetry, etc.) behind these particles. In contrast. 
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collider experiments may identify a long-lived, weakly interacting particle, 

but will not be able to test its cosmological stability or abundance. Only 
by combining the information provided by many different experimental ap- 
proaches is the mystery of dark matter's particle nature likely to be solved. 
Although a confirmed detection of dark matter in any one search channel 
would constitute a discovery of the utmost importance, it would almost 
certainly leave many important questions unanswered as well. 
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